A genotyping array for the globally invasive vector mosquito, Aedes albopictus.

BACKGROUND: Although whole-genome sequencing (WGS) is the preferred genotyping method for most genomic analyses, limitations are often experienced when studying genomes characterized by a high percentage of repetitive elements, high linkage, and recombination deserts. The Asian tiger mosquito (Aedes albopictus), for example, has a genome comprising up to 72% repetitive elements, and therefore we set out to develop a single-nucleotide polymorphism (SNP) chip to be more cost-effective. Aedes albopictus is an invasive species originating from Southeast Asia that has recently spread around the world and is a vector for many human diseases. Developing an accessible genotyping platform is essential in advancing biological control methods and understanding the population dynamics of this pest species, with significant implications for public health. METHODS: We designed a SNP chip for Ae. albopictus (Aealbo chip) based on approximately 2.7 million SNPs identified using WGS data from 819 worldwide samples. We validated the chip using laboratory single-pair crosses, comparing technical replicates, and comparing genotypes of samples genotyped by WGS and the SNP chip. We then used the chip for a population genomic analysis of 237 samples from 28 sites in the native range to evaluate its usefulness in describing patterns of genomic variation and tracing the origins of invasions. RESULTS: Probes on the Aealbo chip targeted 175,396 SNPs in coding and non-coding regions across all three chromosomes, with a density of 102 SNPs per 1 Mb window, and at least one SNP in each of the 17,461 protein-coding genes. Overall, 70% of the probes captured the genetic variation. Segregation analysis found that 98% of the SNPs followed expectations of single-copy Mendelian genes. Comparisons with WGS indicated that sites with genotype disagreements were mostly heterozygotes at loci with WGS read depth < 20, while there was near complete agreement with WGS read depths > 20, indicating that the chip more accurately detects heterozygotes than low-coverage WGS. Sample sizes did not affect the accuracy of the SNP chip genotype calls. Ancestry analyses identified four to five genetic clusters in the native range with various levels of admixture. CONCLUSIONS: The Aealbo chip is highly accurate, is concordant with genotypes from WGS with high sequence coverage, and may be more accurate than low-coverage WGS.

Longitudinal Study of Viral Diversity Associated with Mosquito Species Circulating in Cambodia.

Arthropod-borne viruses (arboviruses) pose a significant global health threat and are primarily transmitted by mosquitoes. In Cambodia, there are currently 290 recorded mosquito species, with at least 17 of them considered potential vectors of arboviruses to humans. Effective surveillance of virome profiles in mosquitoes from Cambodia is vital, as it could help prevent and control arbovirus diseases in a country where epidemics occur frequently. The objective of this study was to identify and characterize the viral diversity in mosquitoes collected during a one-year longitudinal study conducted in various habitats across Cambodia. For this purpose, we used a metatranscriptomics approach and detected the presence of chikungunya virus in the collected mosquitoes. Additionally, we identified viruses categorized into 26 taxa, including those known to harbor arboviruses such as Flaviviridae and Orthomyxoviridae, along with a group of viruses not yet taxonomically identified and provisionally named "unclassified viruses". Interestingly, the taxa detected varied in abundance and composition depending on the mosquito genus, with no significant influence of the collection season. Furthermore, most of the identified viruses were either closely related to viruses found exclusively in insects or represented new viruses belonging to the Rhabdoviridae and Birnaviridae families. The transmission capabilities of these novel viruses to vertebrates remain unknown.

Phylogeny and spatial distribution of Japanese encephalitis virus vector species in Cambodia.

In Southeast Asia, despite the use of Japanese encephalitis vaccines and vaccination coverage, Japanese encephalitis (JE) transmission is still a major public health issue. The main vectors of this virus are mosquitoes from the genus Culex, which diversity and density are important in Southeast Asia. The main vector species of Japanese encephalitis virus (JEV) in Cambodia belong to the Vishnui subgroup. However, their morphological identification solely based on the adult stage remains challenging, making their segregation and detection difficult. In order to identify and describe the distribution of the three main JEV vector species in Cambodia, namely Culex vishnui, Cx. pseudovishnui and Cx. tritaeniorhynchus, mosquito samplings were carried out throughout the country in different environments. Phylogenetic analysis of the cytochrome c oxidase subunit I (coI) gene using maximum-likelihood tree with ultrafast bootstrap and phylogeographic analysis were performed. The three main Culex species are phylogenetically separated, and represent two distinct clades, one with Cx. tritaeniorhynchus and the second with Cx. vishnui and Cx. pseudovishnui, the latter appearing as a subgroup of Cx. vishnui. The phylogeographic analysis shows a distribution of the Vishnui subgroup on the entire Cambodian territory with an overlapped distribution areas leading to a sympatric distribution of these species. The three JEV vector species are geographically well-defined with a strong presence of Cx. pseudovishnui in the forest. Combined with the presence of Cx. tritaeniorhynchus and Cx. vishnui in rural, peri-urban, and urban areas, the presence of JEV-competent vectors is widespread in Cambodia.

A genetic switch controls Pseudomonas aeruginosa surface colonization.

Efficient colonization of mucosal surfaces is essential for opportunistic pathogens like Pseudomonas aeruginosa, but how bacteria collectively and individually adapt to optimize adherence, virulence and dispersal is largely unclear. Here we identified a stochastic genetic switch, hecR-hecE, which is expressed bimodally and generates functionally distinct bacterial subpopulations to balance P. aeruginosa growth and dispersal on surfaces. HecE inhibits the phosphodiesterase BifA and stimulates the diguanylate cyclase WspR to increase c-di-GMP second messenger levels and promote surface colonization in a subpopulation of cells; low-level HecE-expressing cells disperse. The fraction of HecE+ cells is tuned by different stress factors and determines the balance between biofilm formation and long-range cell dispersal of surface-grown communities. We also demonstrate that the HecE pathway represents a druggable target to effectively counter P. aeruginosa surface colonization. Exposing such binary states opens up new ways to control mucosal infections by a major human pathogen.

pWCP is a widely distributed and highly conserved Wolbachia plasmid in Culex pipiens and Culex quinquefasciatus mosquitoes worldwide.

Mosquitoes represent the most important pathogen vectors and are responsible for the spread of a wide variety of poorly treatable diseases. Wolbachia are obligate intracellular bacteria that are widely distributed among arthropods and collectively represents one of the most promising solutions for vector control. In particular, Wolbachia has been shown to limit the transmission of pathogens, and to dramatically affect the reproductive behavior of their host through its phage WO. While much research has focused on deciphering and exploring the biocontrol applications of these WO-related phenotypes, the extent and potential impact of the Wolbachia mobilome remain poorly appreciated. Notably, several Wolbachia plasmids, carrying WO-like genes and Insertion Sequences (IS), thus possibly interrelated to other genetic units of the endosymbiont, have been recently discovered. Here we investigated the diversity and biogeography of the first described plasmid of Wolbachia in Culex pipiens (pWCP) in several islands and continental countries around the world-including Cambodia, Guadeloupe, Martinique, Thailand, and Mexico-together with mosquito strains from colonies that evolved for 2 to 30 years in the laboratory. We used PCR and qPCR to determine the presence and copy number of pWCP in individual mosquitoes, and highly accurate Sanger sequencing to evaluate potential variations. Together with earlier observation, our results show that pWCP is omnipresent and strikingly conserved among Wolbachia populations within mosquitoes from distant geographies and environmental conditions. These data suggest a critical role for the plasmid in Wolbachia ecology and evolution, and the potential of a great tool for further genetic dissection and possible manipulation of this endosymbiont.

Mosquito diversity (Diptera: Culicidae) and medical importance in four Cambodian forests.

BACKGROUND: A total of 290 mosquito species are recorded in Cambodia among which 43 are known vectors of pathogens. As Cambodia is heavily affected by deforestation, a potential change in the dynamic of vector-borne diseases (VDBs) could occur through alteration of the diversity and density of sylvatic vector mosquitoes and induce an increase in their interactions with humans. Understanding mosquito diversity is therefore critical, providing valuable data for risk assessments concerning the (re)emergence of local VBDs. Consequently, this study mainly aimed to understand the spatial and temporal distribution of sylvatic mosquito populations of Cambodia by determining which factors impact on their relative abundance and presence. METHODS: A study was conducted in 12 sites from four forests in Cambodia. All mosquitoes, collected during the dry and rainy seasons, were morphologically identified. The diversity and relative density of mosquito species in each site were calculated along with the influence of meteorological and geographical factors using a quasi-Poisson generalized linear model. RESULTS: A total of 9392 mosquitoes were collected belonging to 13 genera and 85 species. The most represented genera were Culex, accounting for 46% of collected mosquitoes, and Aedes (42%). Besides being the most abundant species, Culex pseudovishnui and Aedes albopictus, which are known vectors of numerous arboviruses, were present in all sites during both dry and rainy seasons. The presence of mosquito species reported to be zoo-anthropophilic feeders was also observed in both forested and urban areas. Finally, this study demonstrated that altitude, temperature and precipitation impacted the abundance of mosquitoes but also influenced species community composition. CONCLUSION: The results indicate an important diversity of mosquitoes in the four forests and an influence of meteorological and geographical factors on their community. Additionally, this work highlights in parallel the abundance of species considered to be of medical importance and therefore underlines the high risk of pathogen emergence/re-emergence in the region.

Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification.

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito rRNA datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNA-seq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Predicting the Effects of Climate Change on Dengue Vector Densities in Southeast Asia through Process-Based Modeling.

BACKGROUND: Aedes aegypti and Ae. albopictus mosquitoes are major vectors for several human diseases of global importance, such as dengue and yellow fever. Their life cycles and hosted arboviruses are climate sensitive and thus expected to be impacted by climate change. Most studies investigating climate change impacts on Aedes at global or continental scales focused on their future global distribution changes, whereas a single study focused on its effects on Ae. aegypti densities regionally. OBJECTIVES: A process-based approach was used to model densities of Ae. aegypti and Ae. albopictus and their potential evolution with climate change using a panel of nine CMIP6 climate models and climate scenarios ranging from strong to low mitigation measures at the Southeast Asian scale and for the next 80 y. METHODS: The process-based model described, through a system of ordinary differential equations, the variations of mosquito densities in 10 compartments, corresponding to 10 different stages of mosquito life cycle, in response to temperature and precipitation variations. Local field data were used to validate model outputs. RESULTS: We show that both species densities will globally increase due to future temperature increases. In Southeast Asia by the end of the century, Ae. aegypti densities are expected to increase from 25% with climate mitigation measures to 46% without; Ae. albopictus densities are expected to increase from 13%-21%, respectively. However, we find spatially contrasted responses at the seasonal scales with a significant decrease in Ae. albopictus densities in lowlands during summer in the future. DISCUSSION: These results contrast with previous results, which brings new insight on the future impacts of climate change on Aedes densities. Major sources of uncertainties, such as mosquito model parametrization and climate model uncertainties, were addressed to explore the limits of such modeling. https://doi.org/10.1289/EHP11068.

dynamAedes: a unified modelling framework for invasive Aedes mosquitoes.

Mosquito species belonging to the genus Aedes have attracted the interest of scientists and public health officers because of their capacity to transmit viruses that affect humans. Some of these species were brought outside their native range by means of trade and tourism and then colonised new regions thanks to a unique combination of eco-physiological traits. Considering mosquito physiological and behavioural traits to understand and predict their population dynamics is thus a crucial step in developing strategies to mitigate the local densities of invasive Aedes populations. Here, we synthesised the life cycle of four invasive Aedes species (Ae. aegypti, Ae. albopictus, Ae. japonicus and Ae. koreicus) in a single multi-scale stochastic modelling framework which we coded in the R package dynamAedes. We designed a stage-based and time-discrete stochastic model driven by temperature, photo-period and inter-specific larval competition that can be applied to three different spatial scales: punctual, local and regional. These spatial scales consider different degrees of spatial complexity and data availability by accounting for both active and passive dispersal of mosquito species as well as for the heterogeneity of the input temperature data. Our overarching aim was to provide a flexible, open-source and user-friendly tool rooted in the most updated knowledge on the species' biology which could be applied to the management of invasive Aedes populations as well as to more theoretical ecological inquiries.

MALDI-TOF MS: An effective tool for a global surveillance of dengue vector species.

Dengue, Zika and chikungunya viruses cause significant human public health burdens in the world. These arboviruses are transmitted by vector mosquito species notably Aedes aegypti and Aedes albopictus. In the Pacific region, more vector species of arboviruses belonging to the Scutellaris Group are present. Due to the expansion of human travel and international trade, the threat of their dispersal in other world regions is on the rise. Strengthening of entomological surveillance ensuring rapid detection of introduced vector species is therefore required in order to avoid their establishment and the risk of arbovirus outbreaks. This surveillance relies on accurate species identification. The aim of this study was to assess the use of the Matrix-Assisted Laser Desorption Ionization Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) as a tool for an international identification and surveillance of these mosquito vectors of arboviruses. Field-mosquitoes belonging to 8 species (Ae. aegypti, Ae. albopictus, Aedes polynesiensis, Aedes scutellaris, Aedes pseudoscutellaris, Aedes malayensis, Aedes futunae and Culex quinquefasciatus) from 6 countries in the Pacific, Asian and Madagascar, were included in this study. Analysis provided evidence that a MALDI-TOF database created using mosquitoes from the Pacific region allowed suitable identification of mosquito species from the other regions. This technic was as efficient as the DNA sequencing method in identifying mosquito species. Indeed, with the exception of two Ae. pseudoscutellaris, an exact species identification was obtained for all individual mosquitoes. These findings highlight that the MALDI-TOF MS is a promising tool that could be used for a global comprehensive arbovirus vector surveillance.

Zika vector competence data reveals risks of outbreaks: the contribution of the European ZIKAlliance project.

First identified in 1947, Zika virus took roughly 70 years to cause a pandemic unusually associated with virus-induced brain damage in newborns. Zika virus is transmitted by mosquitoes, mainly Aedes aegypti, and secondarily, Aedes albopictus, both colonizing a large strip encompassing tropical and temperate regions. As part of the international project ZIKAlliance initiated in 2016, 50 mosquito populations from six species collected in 12 countries were experimentally infected with different Zika viruses. Here, we show that Ae. aegypti is mainly responsible for Zika virus transmission having the highest susceptibility to viral infections. Other species play a secondary role in transmission while Culex mosquitoes are largely non-susceptible. Zika strain is expected to significantly modulate transmission efficiency with African strains being more likely to cause an outbreak. As the distribution of Ae. aegypti will doubtless expand with climate change and without new marketed vaccines, all the ingredients are in place to relive a new pandemic of Zika.

Small-scale field assessment against the dengue vector Aedes aegypti using the auto-dissemination approach in an urban area of Vientiane, Lao PDR.

BACKGROUND: In Lao PDR, dengue fever is the most important vector borne disease and vector control remains the principal method to fight against Aedes aegypti the primary transmitter mosquito species. Vector control management programs need new strategies in addition to conventional larviciding and adulticiding interventions in the country. In this study, we examined the In2Care® Mosquito Trap's efficacy using insecticide auto-dissemination strategy. The insecticide pyriproxyfen, present in powder form inside the trap station, contaminates the body of gravid female mosquitoes visiting the traps and is later on disseminated via the mosquitoes in breeding sites surrounding the traps. We tested the attractiveness of the Traps, their efficacy to reduce the larval and adult abundance, and the impact on emergence rates. Specifically, we tested if the servicing interval of the In2Care® Mosquito Trap could be extended to 12 weeks. METHODS: Two black plastic ovitrap buckets and two BG® sentinel traps were placed in the premises of the Science campus of Vientiane Capital located in an urban area to measure weekly the larval and adult relative abundance of Aedes mosquitoes from 2017 to 2019. Twenty-five In2Care® Mosquito Traps were evenly distributed in this area and two studies of 12 weeks were implemented during January and April 2018 and, July to October 2018 (dry and rainy season, respectively). Every 2 weeks, water samples from 5 In2Care® Traps were randomly selected and tested at the laboratory with Ae. aegypti larvae to measure the larval and pupal mortality. The relative abundance of Aedes mosquitoes in the BG traps® with the presence of In2Care® Traps in 2018, was compared with the surveillance results obtained in 2017 and 2019 without In2Care® Traps. Every week, water samples from the ovitrap buckets were tested for Emergence Inhibition (EI). RESULTS: The In2Care® Traps were very attractive to gravid Ae. aegypti mosquitoes specifically during the rainy seasons with 96% of the traps colonized with larvae/pupae within four weeks. The bioassays showed 100% mortality in the water samples from the traps during the twelve weeks studies showing the good efficacy over time of the pyriproxyfen without additional servicing in the 12 week period. In addition, the larvicide was successfully disseminated into the ovitrap buckets placed in the treated area where 100% of EI during all weeks of intervention was measured. There was no significant effect of the treatment on adult abundance reduction in the treated area, probably due to recolonization of adult mosquitoes surrounding the field experiment. CONCLUSIONS: The observed potential of the In2Care® Mosquito Trap using the auto-dissemination strategy could lead to the use of this new tool in combination with conventional control methods against Dengue vectors in urban tropical areas. Large scale field trials should be implemented in Lao PDR to prove its efficacy for Public Health programs.

Culex vishnui (Diptera: Culicidae): An Overlooked Vector of Arboviruses in South-East Asia.

Culex vishnui Theobald, 1901, a main vector of Japanese encephalitis virus (JEV), is widely distributed in the Oriental region where it often accounts for a great part of the culicid fauna. This species also has been found naturally infected with at least 13 other arboviruses of medical and veterinary importance. Females blood feed predominantly upon pigs and birds, but may readily bite cattle and humans. Because of its abundance, medical importance, and presence throughout ecological gradients among urban, peri-urban, and rural areas, Cx. vishnui potentially may serve as a bridge vector transmitting viruses from natural and wild hosts to humans. Being zoo- and anthropophagic, omnipresent in the Oriental region, and presenting strong resistance to many insecticide families, this overlooked mosquito species may pose a serious health risk in one of the most densely populated regions of the world.

Genomic Shifts, Phenotypic Clines, and Fitness Costs Associated With Cold Tolerance in the Asian Tiger Mosquito.

Climatic variation is a key driver of genetic differentiation and phenotypic traits evolution, and local adaptation to temperature is expected in widespread species. We investigated phenotypic and genomic changes in the native range of the Asian tiger mosquito, Aedes albopictus. We first refine the phylogeographic structure based on genome-wide regions (1,901 double-digest restriction-site associated DNA single nucleotide polymophisms [ddRAD SNPs]) from 41 populations. We then explore the patterns of cold adaptation using phenotypic traits measured in common garden (wing size and cold tolerance) and genotype-temperature associations at targeted candidate regions (51,706 exon-capture SNPs) from nine populations. We confirm the existence of three evolutionary lineages including clades A (Malaysia, Thailand, Cambodia, and Laos), B (China and Okinawa), and C (South Korea and Japan). We identified temperature-associated differentiation in 15 out of 221 candidate regions but none in ddRAD regions, supporting the role of directional selection in detected genes. These include genes involved in lipid metabolism and a circadian clock gene. Most outlier SNPs are differently fixed between clades A and C, whereas clade B has an intermediate pattern. Females are larger at higher latitude yet produce no more eggs, which might favor the storage of energetic reserves in colder climate. Nondiapausing eggs from temperate populations survive better to cold exposure than those from tropical populations, suggesting they are protected from freezing damages but this cold tolerance has a fitness cost in terms of egg viability. Altogether, our results provide strong evidence for the thermal adaptation of A. albopictus across its wide temperature range.

Monitoring insecticide resistance of adult and larval Aedes aegypti (Diptera: Culicidae) in Phnom Penh, Cambodia.

BACKGROUND: Dengue fever is a major public health concern in Cambodia, with thousands of cases every year in urban, suburban and rural areas of the country. The main vector of dengue fever in Cambodia is Aedes aegypti. The organophosphate larvicide temephos and adulticides belonging to the pyrethroid family have been widely used for decades by public health authorities to fight dengue vectors, but resistance of Ae. aegypti to these insecticides has been previously described for Cambodia. METHODS: In order to adapt the vector control strategy presently used in Cambodia, we tested 14 adulticides belonging to the carbamate, organochlorine, organophosphate, and pyrethroid insecticide families and three larvicides [temephos, spinosad and Bacillus thuringiensis ser. israelensis (Bti)] belonging to three different insecticide families (organophosphates, spinosyns and entomopathogenic bacteria). The standard procedures used here to test the adults and larvae of an Ae. aegypti population from Phnom Penh followed World Health Organization guidelines. RESULTS: For adults, high mortality rates were observed with carbamate, organophosphate and organochlorine (with the exception of dichlorodiphenyltrichloroethane) insecticides (i.e. between 87.6 and 100%), while low mortality rates were observed with all of the tested pyrethroid insecticides (i.e. between 1 and 35%). For larvae, no resistance against Bti was detected [resistance ratio (RR90 < 1.6)], but moderate resistance was observed for temephos and spinosad (RR90 < 5.6). CONCLUSIONS: The results of this study indicate that (i) Bti should be considered a serious alternative to temephos for the control of Ae. aegypti larvae; and (ii) the carbamate adulticides propoxur and bendiocarb should be employed instead of the widely used pyrethroid insecticides for the control of adult Ae. aegypti on land under mosaic farming and crop rotation in Cambodia, as the insects were found to be resistant to the latter types of insecticide. Research focusing on insecticide resistance and innovative and effective vector control strategies should be undertaken as a public health priority in Cambodia.

Mosquito (Diptera: Culicidae) diversity and medical importance in Koh Kong mangrove forests, Cambodia.

Background: Mangroves are an ecosystem interface between land and sea, forming distinctive shallow-water marine communities in tropical and subtropical waters. The mangrove forest surface in Cambodia is being reduced due to deforestation. Because the mangrove type of ecosystem generally hosts a great diversity of mosquitoes, the urbanization of these ecosystems will increase interactions between humans and wild mosquitoes, and might thus serve as a potential source of new infectious diseases. Understanding mosquito diversity and analyzing their virome is critical to estimate the risk of emergence or future outbreaks of mosquito-borne diseases. Objective: To understand the mosquito diversity of mangrove forests of Koh Kong province (Cambodia). Methods: In 2019, the mosquito fauna was sampled for 3 consecutive days using BG-Sentinel and light traps, in 3 locations in the mangrove forests of Koh Kong province (Cambodia) during both dry and rainy seasons. Results: A total of 3107 samples were collected, belonging to 10 genera for 34 species. The Culex genus was the most diverse, accounting for 10 species. One species, Culex sitiens, represented over 60% of all collected mosquitoes. A total of 12 medically important species were recorded, 2 species, Aedes (Stegomyia) albopictus and Culex vishnui, were collected in all sites and during both the dry and rainy seasons, highlighting a potential risk of these species acting as bridge vectors. Conclusions: If new arboviruses were to be recorded in this peculiar area, it would indicate that the mosquito species found have the potential to act as a bridge between sylvatic and anthropogenic arboviruses.

High Rickettsial Diversity in Rodents and Their Ectoparasites From the Central Highlands of Madagascar.

Rickettsioses are among emerging infectious diseases around the world. In Madagascar, little information is available regarding Rickettsia (Rickettsiales: Rickettsiaceae) diversity and their potential impacts on public health. In fact, molecular screening of ectoparasites of mammals reported the presence of three species, Rickettsia africae, Rickettsia typhi, and Rickettsia felis. The present study aims to investigate the diversity of Rickettsia in small mammals and associated ectoparasites (fleas and ticks) using a molecular approach. In September and December 2016, fieldworks were undertaken in two districts of Madagascar to capture small mammals using standard traps (Tomahawk and Sherman traps) and collect associated ectoparasites. In total, 12 taxa of ectoparasites (5 flea and 7 tick species) were collected from 89 individuals of four species of terrestrial small mammals. Rickettsia spp. were molecularly identified in one specimen of Rattus rattus (Rodentia: Muridae), one specimen of Pulex irritans (Siphonaptera: Pulicidae) as well as four specimens of Ixodes cf. colasbelcouri (Ixodida: Ixodidae). This study showed the presence of three phylogenetically distinct taxa of Rickettsia in small mammals and their ectoparasites. The current study broadens our knowledge on the diversity of Rickettsia in the Central Highlands of Madagascar and highlights for the first time the presence of Ri. felis in R. rattus and in tick, I. cf. colasbelcouri in Madagascar. Additional studies are needed to have exhaustive information on Rickettsia in small mammals and their ectoparasites, to determine their pathogenicity as well as their potential effects on public health in order to update the national policy for the control of emerging infectious diseases in Madagascar.